Bistable optical response of a nanoparticle heterodimer: Mechanism, phase diagram, 

and switching time 
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We conduct a theoretical study of the bistable optical response of a nanoparticle heterodimer 
comprised of a closely spaced semiconductor quantum dot and metal nanoparticle. The bistable 
nature of the response results from the interplay between the quantum dot's optical nonlinearity 
and its self-action (feedback) originating from the presence of the metal nanoparticle. We show that 
the feedback is governed by a complex valued coupling parameter G. Both the real and imaginary 
parts of G (Gr and G\) play an important role in the occurrence of bistability, which is manifested 
in an S-shaped dependence of the quantum dot excited state population on the intensity of the 
external field, and hysteresis of the population. From our calculations, we find that at Gr = 0, the 
critical value for bistability to occur is Gi = 8r, whereas at Gi = 0, the critical value of Gr = 4F, 
where T is the polarization dephasing rate. Thus, there exist two different (limiting) mechanisms 
of bistability, depending on whether Gr is much larger or much smaller than Gi. We also calculate 
the bistability phase diagram within the system's parameter space: spanned by Gr, Gi and A, the 
latter being the detuning between the driving frequency and the transition frequency of the quantum 
dot. Additionally, switching times from the lower stable branch to the upper one (and vise versa) 
are calculated as a function of the intensity of the driving field. We show that the conditions for 
bistability to occur can be realized, for example, for a heterodimer comprised of a closely spaced 
CdSe (or CdSe/ZnSe) quantum dot and a gold nanosphere. 

PACS numbers: 78.67.-n 73.20.Mf 85.35.-p 



I. INTRODUCTION 

Optical bistability is a fascinating nonlinear phe- 
nomenon, the essence of which is controlling the flow of 
light by light itself. It is of great importance for opti- 
cal technologies, in particular, for optical logic and sig- 
nal processing. The key ingredients for bistable response 
to occur are optical nonlinearity of the material and a 
positive feedback. Interplay of the two can result in a 
multi-valued nonlinear output within a certain range of 
the system parameter space. A generic optical bistable 
clement exhibits two stationary stable states for the same 
input intensity, a property which, in principle, opens the 
door to applications such as all-optical switches, optical 
transistors, and optical memories. 

The phenomenon of optical bistability was predicted 
by McCalli in 1974 and demonstrated experimentally for 
the first time in 1976 by Gibbs, McCall, and Venkatesar* 2 " 
(see also Refs. for an overview) . A Fabry- Perot cavity 
with potassium atoms was used to verify the effect^ It 
has been demonstrated that cavities filled with semicon- 
ductor materials as well as semiconductor micro cavities 
can reveal similar behavior^— 

A vast amount of literature has been devoted to ex- 
plore the topic (an extended bibliography can be found 
in Ref . [Toh . especially on the micro- and nanoscale. Re- 
cently, it was suggested that heterodimers of a closely 
spaced semiconductor quantum dot (SQD) and metal 



nanoparticle (MNP) would be interesting nanoscale sys- 
tems that exhibit bistable optical response! 11 ' 12 In fact, 
such systems have a variety of interesting optical prop- 
erties. In particular, control of the SQD exciton emis- 
sio n 13 i 14 and an increased exciton-hole recombination 
ratep^ variable quenching of the SQD photoluminescence 
by proximate gold nanoparticles^ nonlinear Fano reso- 
nances ; 11 ' 17 optical hysteresis of the Rayligh scattering^ 
"meta-molecular" resonances^ inhibition of optical ex- 
citation and enhancement of Rabi flopping^ tunable 
nanoswitching^i and gain without inversion 2 ^ have been 
predicted and partly observed^ The role of the multi- 
pole SQD-MNP interaction in explaining the spectra of 
hybrid systems has been discussed in Ref. H^. All these 
effects are governed by the geometrical and material pa- 
rameters of the hybrid cluster, thus providing a promising 
opportunity of controlling spectral and dynamical behav- 
ior of nanoscale objects. The optical response of SQD- 
MNP heteroclusters result from strong coupling between 
excitons in the SQD and plasmons in the MNP. This 
coupling may be tuned by varying the system's materials 
and geometry, giving rise to a dramatic modification of 
the response. The novel properties exhibited by such hy- 
brid nanosystems hold the potential to revolutionaze op- 
toelectronics and nanophotoniesi 24 ' 25 Additionally, new 
materials (frequently called metamaterials) such as pho- 
tonic crystals^ surface-plasmon polaritonic crystals^ 
and materials with a negative index of refraction^ have 
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been shown to manifest bistable behavior. 

In this paper, we present an important step in a fur- 
ther understanding of the optical response of an SQD- 
MNP hcterodimcr. We analyze the problem in all its 
facets, adding to previous work by other author s 11 ' 12 ' 18 
a comprehensive analysis of the system's parameter sub- 
space where bistability may occur (the so-called phase 
diagram), examples of realistic conditions under which 
bistability may actually be achieved with existing mate- 
rials (CdSc quantum dot and gold nanoparticlc at various 
distances), a fundamental understanding of the mecha- 
nism of bistability, and a study of the switching of the 
system between both stable branches. 

With regards to the mechanism, we focus on the role 
of the SQD-MNP (complex) coupling parameter G = 
Gr + iGj, which quantifies the self-action (feedback) for 
the SQD in the presence of the MNP. We distinguish 
between the roles of the real and imaginary parts of G 
(Gr and Gi) and show that they result in two differ- 
ent mechanisms of the SQD bistability. In the case of 
Gr, 7^ and G\ = 0, the feedback is provided by the 
population-dependent resonance frequency of the SQD, 
while in the other case, G\ ^ and Gr = 0, it originates 
from the destructive interference of the driving field with 
the secondary field produced by the SQD. We calculate 
the bistability phase diagram within the system's param- 
eter space: spanned by Gr, Gi and A, the latter being 
the detuning between the driving frequency and the tran- 
sition frequency of the quantum dot, and uncover a pe- 
culiar behavior of the bistability threshold as a function 
of Gr and G\. The switching time between both stable 
branches is calculated as a function of intensity of the 
driving field, which is important from the viewpoint of 
practical applications as an all-optical switch. 

This paper is organized as follows. In the next section, 
we present the system setup and analyze the fields ex- 
perienced by the SQD and the MNP, both exposed to a 
driving field. Sec. Ill deals with the density matrix for- 
malism for describing the optical dynamics of the SQD 
coupled to the MNP. In Sec. IV, we discuss in detail 
the conditions for bistability to occur in the SQD optical 
response, based on calculations of the bistability phase 
diagrams. The physical interpretation of the influence of 
the SQD-MNP coupling parameter G = Gr + G\ and the 
detuning away from the SQD resonance is presented. In 
Sec. V, we study the switching time of the system when 
subjected to a sudden change in the driving intensity. In 
Sec. VI, we summarize and conclude. 



II. SYSTEM SETUP 

The system of our interest is schematically shown in 
Fig. [TJ It is comprised of a single spherical SQD, char- 
acterized by a bare dielectric constant e s , coupled to 
a closely positioned spherical MNP with polarizability 
a(uj). This heterodimer is assumed to be embedded 
in a dielectric background with dispersionless permit- 



tivity Sb and driven by a monochromatic external field 
E = (1/2)Eq exp(— iuot) + c.c. which is linearly polar- 
ized along the SQD-MNP axis. The frequency of the 
incident field w is assumed to be close to the bare ex- 
citon transition frequency loq which, in turn, is close to 
the plasmon resonance peak wsp. We denote the radii 
of the MNP and the SQD as a and r, respectively, while 
the center-to-center distance between the particles is d. 
These three parameters (a, r, and d) are assumed to be 
small as compared to the SQD emission wavelength, al- 
lowing us to neglect retardation effects and to consider 
both nanoparticles as point dipoles. 




FIG. 1. Schematics of a SQD-MNP heterodimer embedded 
in a homogeneous dielectric host with permittivity Eb and 
subjected to an external field of amplitude Eq, polarized along 
the system axis. e s is the SQD bare dielectric constant, a(ui) 
is the polarizability of the MNP. -Bsm and -Ems, respectively, 
denote the electric fields produced by the polarization of the 
SQD at the position of the MNP and vice versa. 

The dominant optical excitations of the SQD are con- 
fined excitons with a discrete energy spectrum. We re- 
strict ourselves to taking into account only one (lowest) 
exciton energy level characterized by a narrow absorption 
line width and a transition dipolc moment fi. The opti- 
cal dynamics of the exciton transition will be described 
quantum mechanically by making use of the Maxwell- 
Bloch equations for the 2x2 density matrix p mn , where 
m and n may be (for the ground state) or 1 (for the 
excited state). 

The MNP is considered classically in the qua- 
sistatic approximation; its response is described by the 
frequency-dependent polarizability a(co) within the point 
dipolc approximation (this can be easily generalized to 
the case of more complex MNP shapes by using an appro- 
priate polarizability tensor). The SQD-MNP interaction 
will be treated within the point dipole-dipole approxima- 
tion. 

Now, let us calculate the fields experienced by the SQD 
and MNP. The external field polarizes the nanoparticles. 
The polarization of the SQD generates an additional field 
Esm at the position of the MNP and vice versa Ems i see 
Fig. [TJ These fields are superposed on the external field 
Eq; so that the fields acting upon the SQD and MNP 



3 



are Eq + Ems and Eq + Eqm, respectively (acting along 
the system axes). Note that all above relationships are 
written for the field amplitudes; the oscillations with the 
optical frequency u> already have been extracted. 

Considering the SQD's induced dipole moment -Psqd 
as a point dipole, the field Esm can be written in the 
form (see, e.g., Refs. |29| andl30h: 
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Here it is assumed that the SQD is a uniformly po- 
larized sphere, the field of which is screened only by the 
dielectric constant e b of the host medium^ Within the 
density matrix formalism, -Psqd = —ifJ>R, where R is 
the amplitude of the off-diagonal density matrix element 
pio = — (i/2)i?exp(— iuit). The MNP dipole moment is 
now determined by the total field E + E$m, i-e, 



-Pvinp = e e b a(uj) Eq + 



a(ui) = Altar 



-PSQD 
27T£ e& d 3 



e m {uS) + 2e b 



(2a) 



(2b) 



Here, e m is the permittivity of the MNP. The peak of the 
MNP polarizability a(u>), when the denominator is min- 
imal, determines the MNP (surface) plasmon resonance. 
We do not take into account the corrections to a due to 
the depolarization shift and radiative damping^ which 
are both negligible for the MNP sizes of our interest (< 
10 ran), 

The field produced by the MNP at the SQD, E U s, 
takes the same form as Eq ([1]), with Psqd replaced by 
PWp- The total field experienced by the SQD equals 
E a + Ems- However, the field inside the SQD should 
be reduced by an effective SQD dielectric constant e' s = 
(e s + 2g;,V/(3efr) (see, e.g., Ref. HH, chapter V, page 138, 
and Ref. [30h . Taking all this into account, the total field 
inside the SQD reads: 
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Eq. (|3|) shows two effects for the SQD due to the presence 
of the MNP. In the first term, one can see a renormal- 
ization of the external field amplitude Eq by a factor 
(l/e' s )[l + a(ui)/ (4ir d 3 )]. The second term reveals a self- 
action of the SQD via the MNP: the field that the SQD 
experiences, depends on its own state through its dipole 
moment amplitude -Psqd- As we will show below, this 
drastically affects the SQD-MNP heterodimer optical re- 
sponse. 

Here, a comment on the second term in Eq. (j3|) is in 
order. In a number of recent publications, dealing with 
the same system, a different formula for this term was 
used, in which the factor e' s in the denominator appeared 
squared . 11 ! 12 ' 14 ^ 7 ' 19 " — We do not agree with this and 
follow the arguments of Ref. [TH that the above factor 
should be linear. 



III. DESCRIBING THE SQD OPTICAL 
DYNAMICS 



As we already mentioned in the previous section, the 
SQD is assumed to be a two level system, having its filled 
valence band as ground state |0) and the lowest exciton 
level as its excited state |1); both states are separated 
by the transition frequency uiq. This approximation is 
justified when the frequency of the external field uj is 
close to the exciton resonance (uj ~ u>o). Throughout 
this paper, we use the rotating-wave approximation, so 
that the time-dependent quantities arc the amplitudes of 
the density matrix elements. The corresponding set of 
equations reads 



z = -i(z 



i) - \ m* 



n*R) 



R = - (T + iA) R + m 



(4a) 



(4b) 



where Z = pn — poo is the population difference between 
the excited and ground states of the SQD and R is the 
amplitude of the off-diagonal density matrix element de- 
fined as p\Q = —(i/2)Rexp(—iujt). The population differ- 
ence Z and the amplitude R of pio are quantities slowly 
varying on the scale of an optical period. The constants 
7 and T represent the rates of population and phase re- 
laxation, respectively, A = uj q — ui is the detuning away 
from the SQD resonance, and f2 = pEsQu/h is the total 
electric field inside the SQD (in frequency units). 

According to Eq. ([3]), the total field ft acting inside the 
SQD can be written in the form 



O = n -iGR 



where f2o and G are given by 



«n = 



G 



1 



a(ui) 
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(5) 



(6a) 



(6b) 



Here, fin is the Rabi frequency of the external field 
renormalized because of the SQD-MNP coupling, with 
On = fJ-Eo/H being the bare Rabi frequency. As we al- 
ready mentioned in the previous section, the second term 
in Eq. ([5]) describes the self-action of the SQD via the 
MNP. The complex-valued constant G = Gr + iG\ is a 
feedback parameter which is determined by the dimer's 
geometry and material properties. Its real part describes 
the near-zone feedback field, while the imaginary part 
is a radiation (far-zone) feedback field (see below). The 
parameter G contains all information governing the SQD 
self-action, such as material constants, geometry of the 
system, and/or details of the interaction (e.g. contribu- 
tions of higher multipoles^) . 



4 



In order to shed light on the effect of self-action on 
the SQD optical dynamics, we substitute Eq. ((3J) into 
Eq. (|4b| and obtain 

R = -[(T-G I Z) + i(A + G R Z)]R + n Z . (7) 

From Eq. (0, it becomes apparent that the SQD self- 
action has two consequences: (i) - the renormalization 
of the SQD resonance frequency ujq — > luq + Gr Z and 
(ii) - the renormalization of the dipole dephasing rate 
r — > r — G\ Z\ both rcnormalizations depend on the pop- 
ulation difference Z . This makes Eqs. (|4a[) and (|4b[) non- 
linear. Similar renormalizations have been reported in re- 
lation with the nonlinear optical response of dense solid 
stat o 32 ' 33 and gaseous 3 ^ assemblies of two-level atoms, 
optically dense thin films , 10 ' 35 " — and linear molecular 
aggregates. 39 ' 40 Both the population dependencies of the 
SQD resonance frequency and the dipole dephasing rate 
provide feedback mechanisms that can give rise to bista- 
bility (see below). 



IV. BISTABILITY OF THE OPTICAL 
RESPONSE 

A. Steady state regime 

First of all, we are interested in steady-state solutions 
of Eqs. (|4"a)l and (|4b|) . which the system reaches after 
turning on the driving held and waiting until the tran- 
sient processes are over. Formally, this can be done by 
setting the time derivatives in Eqs. (|4a| and (|4b|) to zero. 
After simple algebra, we obtain: 

|»o| 2 Z + l Kr-d^ + ^A + GnZ)! 2 

7 r z r 2 ' 



(T-GiZ)+i(A + G R Z) ■ y ' 

As is seen, Eq. (|8aj) is a closed equation which is of 
third order in Z . This implies that, depending on the 
values for A, 7, T, and G it may have three real solu- 
tions. The same applies to the dipole moment amplitude 
R. It should be noticed that the possibility of having a 
three- valued solution to Eq. [5a] implies three- valued op- 
tical response of the SQD-MNP hybrid dimer. 

B. CdSe-Au heterodimer 

Before performing a general study of Eq.( [SaJ) (see 
Sec. IIVC]) . we present the results of its solution obtained 
for a specific system consisting a CdSe SQD coupled to 
a gold MNP. In our numerical calculations, the following 
set of the SQD parameters was selected: the transition 
energy huo = 2.36 eV (which corresponds to the optical 



transition in a 3.3 nm radius SQD), the transition dipole 
moment (i = 0.65 e- nm, the SQD bare dielectric con- 
stant e s = 6.2, the host dielectric permittivity £6 = 1, 
and the SQD relaxation constants 7 = 1.25 ns _1 and 
r = 3.33 ns" 1 ^ We chose the MNP radius a = 10 nm 
and the MNP-SQD center-to-center distance d = 17 nm. 
We have used the tabulated data for the permittivity of 
gold e m (w)^i to calculate the MNP polarizability a(uj). 
according to Eq. (|2b|). We found that a(u>) has a peak 
at Hlusp = 2.4 eV with a width on the order of 0.25 eV 
(see also Ref . [T3) ■ These data allowed us to extract the 
feedback parameter G using Eq. (|6b[) . As is seen from 
that equation, G is a function of frequency. However, 
the frequency domain of our interest is determined by a 
narrow region around the SQD sharp resonance (at most 
of the order of 10 T, see below), whereas the MNP plas- 
mon peak is much broader. Therefore, G is required just 
at the SQD resonance frequency loq. At this frequency, 
for the set of parameters used, G = (27.1 + ll.li)r. 

Fig. [2] shows the results of our calculations together 
with the SQD level diagram hinting at the mechanism of 
the SQD response. Calculations have been performed for 
the set of parameters described above, choosing a bare 
SQD detuning A = IT (i.e. the incident field is almost in 
resonance with the bare SQD transition). However, when 
the incoming field amplitude is small, so that Z « — 1, 
the SQD resonance is red-shifted by the MNP to wo — Gr 
[left panel of Fig.[2][b)]. The renormalizcd detuning A' = 
A — Gr = — 26. IT then is pretty large. In other words, 
for small f^o, the SQD is strongly out of resonance. 

In Fig. [H^a), the colored solid lines (red, black, and 
blue) represent different branches of the S-shapcd char- 
acteristics of Z (population difference) versus |fio| 2 /(7r) 
(dimcnsionless input intensity), clearly indicating the ex- 
istence of a three- valued solution of Eq. (|8a|) . In order to 
determine the stability of different branches, we made use 
of a time domain analysis based on Eqs. (|4"aj) and (|4b|) . 
sweeping adiabatically up and down the external field in- 
tensity |f2o| 2 /(7r). The picture of the system's behavior 
is as follows. On increasing |i7o | 2 / (t^) from zero (indi- 
cated by red arrows), the SQD follows the lower (red) 
branch of the S-shaped curve. Simultaneously, the pop- 
ulation difference Z grows (become less negative), thus 
bringing the SQD closer to resonance with the external 
field [see left panel of Fig. [UJb)]. This holds until the 
intensity reaches a (right) critical magnitude. Upon fur- 
ther increasing |f2o| 2 /(7r), the SQD suddenly jumps to 
the upper (blue) branch of the S-shaped curve, where 
the population difference Z « 0, meaning that the SQD 
is saturated; the heterodimer is now almost in resonance 
with the external field. 

Upon sweeping back the incoming intensity |f2 | 2 /( 7 r) 
[indicated by blue arrows in Fig. [Ha)], the SQD follows 
the upper (blue) branch of the S-shaped curve, keep- 
ing the SQD saturated (Z w 0), until its value reaches 
the other (left) critical magnitude. Further decreasing 
l^o | 2 / (7r) leads to a sudden jump of the SQD population 
difference to the lower (red) branch of the S-shapcd curve, 
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FIG. 2. (a) The steady state solution of Eq. [8a] for a CdSe- 
Au heterodimer, with SQD radius 3.3 nm, MNP radius 10 
nm, and SQD-MNP center-to-center distance 17 nm. Red, 
black, and blue branches form an S-shaped characteristics of 
the SQD population difference Z versus the dimensionless ex- 
ternal field intensity \Q,o\ 2 /yT calculated for a bare detuning 
A = ir. The arrows show the system's time domain route 
(hysteresis loop) obtained after solving Eqs. I4al and I4bl under 
adiabatic sweeping up and down of the external field intensity, 
(b) Schematics of the renormalized SQD energy levels (hor- 
izontal dashed lines) when holding the system at the lower 
(left) and upper (right) stable branches. The green vertical 
arrows indicate the position of the incident frequency with re- 
spect to the SQD energy levels. Red and blue vertical arrows 
show the direction of moving the SQD energy levels when 
sweeping up and down the driving intensity |f2o| 2 /7r. 



thus closing the hysteresis loop. At the lower branch, the 
external field is again far out of resonance with the SQD. 

From these calculations, we conclude that only the 
lower and upper branches of the S-shaped characteris- 
tics are stable. The intermediate branch (black curve 
with a negative gradient) is unreachable (unstable) when 
adiabatically sweeping the incoming field intensity. This 
allows us to speak of bistability in the optical response 
of the SQD-MNP heterodimer. 



C. General study: Bistability phase diagram 

It is of interest to perform a general study of the sys- 
tem's bistability, examining the occurrence of the effect 
in the parameter space Gr, G\, A, and T. As follows 
from Eq. (|8a|) . the relaxation constant T can be used as 
a unit for Gr, Gt, and A and thus is not a relevant pa- 



rameter. 

Our study is based on Eq (|5aj) which, as was already 
mentioned in Sec. IIV Al is of third order in Z. There- 
fore, this equation may have three real roots for specific 
values of Gr, Gt, and A. These three solutions are dif- 
ferent when |f2n| 2 /(7r) in Eq. (|5a|) . formally considered 
as a function of Z, has a minimum and a maximum. 
The threshold for bistability is determined by the con- 
dition that the derivative of |firj| 2 /(7r) with respect to 
Z has a degenerate root (merged extrema). We used 
this definition to calculate the bistability phase diagram. 
The results are shown in Fig. [3J Colored areas denote 




5 10 15 20 
A[r] 





G R [F] 



FIG. 3. Bistability phase diagrams of the SQD optical re- 
sponse. Colored areas mark the parameter subspace where 
bistability exists. From left to right: the phase diagram in the 
parameter subspace (a) [Gr; A; Gi = 0], (b) [Gi; A; Gr = 0], 
and (c) [Gi; Gr]. In the last diagram, the point (Gr; Gi) falls 
in the bistability region of there is at least one value of A for 
which (Gr; Gi; A) yields a bistable response. The boundaries 
between white and colored regions represent the bistability 
threshold within the corresponding parameter sub-space. 

the parameter sub-space where bistability occurs. The 
boundaries between white and colored regions represent 
the bistability threshold for given parameters. One can 
conclude, observing the plots presented in Fig. |3l that 
both Gr and G\ play an important role in the occur- 
rence of bistability. 

Fig. Eta) shows the bistability phase diagram within 
the parameter sub-space [Gr;A;Gt = 0]. First, we ob- 
serve that there is an absolute threshold for bistability 
to occur with respect to Gr: the effect exists only if 
Gr > 4r. This is in agreement with the analytical result 
derived by Fricdberg et al.— for a dense gaseous medium. 
Also, not surprisingly, it is seen that increasing the detun- 
ing A requires a larger Gr to get a bistable response. The 
strong asymmetry of the phase diagram [Gr; A; Gi = 0] 
with respect to a change of sign of A has a natural ex- 
planation. Remember that the mechanism of bistabil- 
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ity in the present case originates from the population- 
dependent SQD resonance frequency shift which brings 
the system into resonance with the incoming field (see se- 
cion lIVB|) . At a positive detuning (cjo > reaching res- 
onance can be achieved within the window < A < Gr: 
changing the population difference — f < Z < allows 
this. At a large negative detuning (uq < lo) bistable be- 
havior cannot be realized because the resonance between 
the the SQD and the external field requires a significant 
positive population difference Z, which is unreachable 
under steady state conditions. 

Within the parameter sub-space [Gi;A;Gr = 0] the 
absolute threshold for bistability is given by Gi = 8r 
[see Fig. [3jb)]. At smaller Gi, the effect is not present 
at all. This result can be verified analytically by set- 
ting Gr = A = in (|5a|) and analyzing the derivative 
of 1 17 1 2 / (t r) , as explained above^ Increasing the de- 
tuning A requires larger values of Gi, however, unlike 
the previous sub-set of parameters [Gr; A; Gi = 0], the 
bistability phase diagram here is symmetric with respect 
to changing the sign of A. The mechanism of bistabil- 
ity in this case is as follows. Using Eq. (|8b|) . it can be 
easily shown, that at a low level of excitation (Z « — I ), 
the feedback field —iGR is out of phase with the ex- 
ternal one f^o, and at G\ 3> T almost compensates the 
latter. The total field f2 = Slo — iGR is on the order of 
r/G i 10 ' 38 , i.e., is very small, thus preventing bistability 
to occur. As the system is being excited, the compen- 
sation decreases, leading to an increase of the total field 
fl inside the SQD and saturating the SQD transition. 
Finally, such an interference-based mechanism also gives 
rise to the second self-sustaining stable state (the upper 
branch), provided G\ is sufficiently large. 

Fig. EJc) shows the bistability phase diagram within 
the parameter sub-space [Gr; Gi; any A] (see caption for 
explanation). The diagram is symmetric with respect 
to the transformation Gr to — Gr. Because of that, we 
present it only for Gr > 0. Remember that the col- 
ored area shows the range of Gr and G\ where bistability 
may exist. Observing the diagram, we can make several 
conclusions. First, the absolute bistability threshold at 
Gi = is Gr = 4r, whereas at Gr = 0, it is Gi = 8T, 
in accordance with the results presented in Figs. Eta) 
and[^b) . Second, within the range < Gr < 4r (be- 
low the bistability threshold with respect to Gr = 4r 
at Gi — 0), the bistability threshold with respect to Gi 
decreases from Gi = 8r to Gi = 6.2T, meaning that 
increasing Gr within this range works towards bistabil- 
ity to occur. Finally, when 4r < Gr < 5.2T (above the 
bistability threshold with respect to Gr = 4r at Gi = 0) , 
there is a range of Gi values, depending on the value of 
Gr, where bistability does not exist. The presence of 
this area in the phase diagram originates from the com- 
plicated interplay of two fields: the renormalized external 
field Slo and the feedback field — iG R, which both deter- 
mine the total field fl inside the SQD, see (O. Within this 
area, these two fields interfere destructively with each 



other, thus preventing the occurrence of bistability. 
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FIG. 4. Steady state solutions of Eq. (|8a[l for various values 
of the center-to-center distance d between the SQD and the 
MNP: d = 40 nm, d = 21.5 nm, d = 18 ran. The set of other 
parameters is the same as in Sec. II V 51 Corresponding values 
of the feedback parameter G are shown in the legend. 

As is seen from Eq. (|6b[) . we can control the strength 
of the feedback parameter G by changing the SQD-MNP 
center-to-center distance d. In Fig. 21 we present the so- 
lutions of Eq. (|5a|) for d = 40, 21.5, and 18 nm. For these 
values of d, the corresponding feedback parameters are 
G = (O.I6 + 0.07t)r, (6.61 + 2. 72i)T, a nd (19 .20 + 7.89 i)T. 
The other parameters are as in Sec. IIV Bl for the CdSe- 
Au dimer. Observing the plots, one can conclude that 
for d = 40 nm [G = (0.16 + 0.07£)r], the dependence 
of the SQD population difference Z on the external field 
intensity |f2o| 2 /(7r) is single- valued (bistability does not 
occur). This is because the magnitudes of Gr and Gi 
are below the bistability threshold [see the phase dia- 
gram in Fig. EI C )]- For the smaller SQD-MNP spacing 
of d = 21.5 nm [G = (6.61 + 2.72i)r], a three-valued so- 
lution of Eq. (|8a[) starts to appear, and the population- 
vs-intensity curve shows an S'-like shape within a narrow 
range for the external field intensity. At this spacing d, 
the feedback parameter magnitudes Gr and Gi are above 
the bistability threshold, as follows from the phase dia- 
gram [Fig. [3]^c)]. A clear S'-like shape of the population- 
vs-intensity dependence is recovered for an even smaller 
SQD-MNP center-to-center distance d = 18 nm; the cor- 
responding G = (19.20 + 7.89i)T is well inside the bista- 
bility region [see again Fig. (3fc)]. 



V. SWITCHING TIME 

In this section, we address the switching time t of the 
SQD-MNP hetcrodimer between the stable branches of 
the bistable Z-vs-Iq characteristics in the vicinity of the 
switching points, i.e. intensities close to the upper and 
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lower critical intensities [marked by red and blue dashes 
lines in Fig.[2Ja)]. Not only it is of fundamental interest 
to investigate the switching dynamics in this nonlinear 
system, it also is of importance in order to assess the 
potential usefulness of such systems as building blocks 
of real devices. Figs. [5] and [6] show the results obtained 
for the upper-driving-intensity switching point [see red 
dashes in Fig. [Ha)]. We defined r as the time which it 
takes for the population difference Z to acquire its first 
maximum after suddenly switching on the incident inten- 
sity, I = |fi | 2 /(7r), to a value slightly larger than the 
critical one I c = \Qo c \ 2 / (l^) . From Fig. [SJ it is clearly 
seen that r sensitively depends on the excess of Jo over 
I c = 531.4125 (the latter is calculated for the set of pa- 
rameters as used in Fig. [2]): the system response drasti- 
cally slows down when the driving intensity Iq approaches 
the critical value I c . Without showing details, we men- 
tion that below the lower critical point [blue dashes in 
Fig. Ufa)], the dynamics of Z approaching its steady- 
state value when starting from Z = and switching off 
the external held is exponential with a characteristic ex- 
ponent on the order of the population relaxation time 

In Fig. IHIa), we plotted the dependence of r (defined 
as explained above) on the excess Jo — Ic of the inci- 
dent intensity Iq over the critical one I c . The numerical 
data points (symbols) can be well fitted by the formula 
r = 5.257 x 10 6 (/ - Ic) -0 ' 505 . It is of great interest 
whether the exponent in the r-vs-J c dependence, approx- 
imately equal to 0.5, is universal. In order to investigate 
this, we performed a series of calculations of the inten- 
sity dependence of the population relaxation time r close 
to the high-intensity switching point, varying the inter- 
particle center-to-center distance d (the coupling param- 
eter G, in other words) and the off-resonance detuning A. 
The results (in log- log scale) are presented in Fig. EJb). 
As is seen from the numerical data and fits, the exponent 

0.5 indeed seems to be universal. The slight deviation 
of the data points from a straight line may be a conse- 
quence of the definition of the relaxation time r as the 
time the population difference Z acquires its first maxi- 
mum after switching on the incident field (see Fig. [5]). 

The kinetics of the population difference Z in the vicin- 
ity to the lower critical point of the Z-vs-Iq bistable char- 
acteristic also show an oscillatory behavior, but not as 
sharp as in the vicinity of the upper critical point. There- 
fore, the definition of the switching time r used above is 
not useful in this case. This can be understood from the 
fact that here the driving intensity is too low to support a 
population difference close to zero. Thus, the population 
relaxation time is dominated by the radiative decay. 



VI. CONCLUSION 

We conducted a theoretical study of the optical re- 
sponse of a nano-heterodimer comprised of a closely 
spaced spherical semiconductor quantum dot and a metal 
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FIG. 5. (a) Time evolution of the SQD population difference 
Z after suddenly switching on the driving field Io above the 
upper critical intensity of the bistable Z-vs-Iq characteristics 
in Fig.Ha), J c = \n 0c \ 2 /(jr) = 531.4125. (b) As (a) but now 
in close proximity to the critical intensity, demonstrating the 
slowing down of the population relaxation when approaching 
the critical intensity from the above [compare the time scales 
in (a) and (b)]. The system parameters were taken as in 
Fig.® 



nanosphere coupled to each other by dipole-dipole forces. 
The coupling results in a self-action of the SQD via 
the MNP, characterized by a complex coupling constant 
G = Gr + iGi, which causes the SQD transition fre- 
quency (through Gr) and dephasing rate (through Gi) 
to depend on the SQD excited state population. This 
provides a feedback mechanism resulting in bistable op- 
tical response of the system (an S-shaped behavior of 
the SQD population difference Z versus incident inten- 
sity I). In particular, at certain field intensities the state 
of the SQD (and the system as a whole) can be switched 
abruptly from a low level of population of the SQD to 
a higher one by changing slightly the input field inten- 
sity. Additionally, the system shows a hysteresis behavior 
when sweeping the incoming intensity up and down. 

The different physical meanings of the coupling con- 
stants Gr and G\ imply two different mechanisms of the 
SQD bistability. If G\ = 0, the feedback is provided 
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FIG. 6. (a) Relaxation time r of the population difference Z 
as a function of the incident intensity To = |fio| 2 /(7r) (above 
the upper critical value I c = |fi c| 2 /(7P) = 531.4125). Pa- 
rameters are the same as as in Fig. The circles are data 
points, while the solid line represents a power-law fit given 



by r = 5.257 x 10 b (T - IcY 



(b) As (a) calculated 



for different center-to-center distances d between the SQD 
and MNP and for different bare detunings A. The other pa- 
rameters are the same as as in Fig. [2] The symbols are the 
data points, while the lines represent the best fits to straight 
lines, given by 0.653 x 10 7 (/ - T c )~ ' 499 (blue line), r = 
5.257x10^ (T -T c )-°' 505 (red line), 0.405xl0- 7 (/ o -T c )-°' 521 



produced by the SQD. Because of that, the thresholds 
for bistability to occur are different in these two cases: 
at G\ = 0, the threshold occurs at Gr = 4r, whereas at 
Gr = 0, it is Gi = 8r. When both constants, Gr and 
Gi, are not zero, the two mechanisms of bistability in- 
terfere, resulting in a quite complicated behavior of the 
bistability threshold as a function of Gr and Gj. We 
calculated the bistability phase diagrams within the sys- 
tem's parameter space: Gr, Gi and A to uncover this 
behavior. Computations performed for a heterodimcr 
comprised of a CdSe quantum dot and an Au nanopar- 
ticle show that bistable behavior may occur for realistic 
nano-particle sizes. It should be noticed that a dimcr 
comprised of strongly coupled two-level molecules can not 
manifest bistability because of the discreteness of the sys- 
tem's levels. The later prevents a continuous change of 
the transition frequency, which is a necessary ingredient 
to achieve positive feedback.^ Thus, a SQD-MNP het- 
erodimcr is a unique nanoscopic system exhibiting this 
feature. Having performed the steady state analysis of 
the optical response of SQD-MNP heterodimers, we also 
studied the time it takes for the system to switch from 
one stable state to the other as a function of the excess of 
the incident intensity with respect to the critical switch- 
ing value. At the upper switching point, we found a 
power-law, which surprisingly, has a universal exponent. 
This switching time diverges when the incident intensity 
Io approaches the critical one I c , indicating the critical 
slowing down of the system response. The switching time 
in the vicinity of the lower critical point does not show 
such pecularities and is on the order of the population 
relaxation time 7 -1 . 

The model we considered is the simplest hybrid nan- 
odimcr. We expect, however, that more complicated 
clusters (such as a SQD surrounded by several MNPs, 
as reported in Ref. [3) or a single quantum dot (or quan- 
tum dot lattice) on top of a metal surface, can also exhibit 
the effects mentioned above; the MNPs just play the role 
of "resonant mirrors" and provide the positive feedback 
which is one of the essential ingredients for bistability to 
occur. 



(black line), and 0.267 x 10~ 7 (T - h 



(green line). 
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